The Size Distribution of Kuiper belt objects for D > 10 km 

Wesley C. Fraser 1,2, JJ. Kavelaars 2,1 

ABSTRACT 

We have performed a survey of the Kuiper belt covering ~ 1/3 a square degree 
of the sky using Suprime-cam on the Subaru telescope, to a limiting magnitude of 
77150 ~ 26.8 and have found 36 new KBOs. We have confirmed that the luminosity 
function of the Kuiper belt must break as previously observed (Bernstein et al. 2004; 
Fuentes & Holman 2008). From the luminosity function, we have inferred the underlying 
size distribution and find that it is consistent with a large object power-law slope qi ~ 
4.8 that breaks to a slope q2 ~ 1.9 at object diameter Db ~ 60 km assuming 6% 
albedos. We have found no conclusive evidence that the size distribution of KBOs with 
inclinations i < 5 is different than that of those with i > 5. We discuss implications of 
this measurement for early accretion in the outer solar system and Neptune migration 
scenarios. 
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Introduction 



The K uiper belt is a pop ulation of planetesimals with diameters as lar ge as a few tho usand 
kilometers (Brown et al. 2006) and a mass of a pproximately a few O.OIM^ ( Gladman et alj [2001: 



Bernstein et al 



2004 



Fuentes HolmanI 120081 ) . The state of the belt is enigmatic, as such large 



objects are not likely to f orm in such a low mass Kuiper belt, over the age of the solar system 
(jKenyon &: Bromley! l200ll ). Rather, it is likely that a much more massive initial belt underwent 
accretion for a t i me before sorne larg e-scale mass-depletion event occurred, such as a stellar passage 



Ida et al 



2000 



Levison et al 



Leyison et alJ 12004 ) . a sweep through of the mea n-motion resonances of Neptune 



200a), or the scattering of KBOs by rogue planets (IGladman fc Chanll2006l ). What 



ever the event, it is also responsible, at least in part, for the dynamically excited orbits of KBO 
populations such as the plutinos or the scattering population ([Gladman et al.ll2008l ). 



The size distribution of the Kuiper belt contains a fossil record of the end-state of the accretion 
processes that occurred in that region. Knowledge of the size distribution can constrain disruption 
strengths of the bodies, formation time-scales in the outer solar system, and the early conditions of 
the proto-planetary disk. This makes the determination of the size distribution a primary constraint 
on Kuiper belt formation scenarios. 

Because of the large distance to the Kuiper belt, the size distribution is not determined directly, 
but rather, is inferred from the shape of the observed luminosity function (LF). Early observations 
determined that the LF for brig ht objects (rn(R) < 26) was well represented by a power-la\y 
with a logarithinic slope a ~ 0.7 djewitt et al.1 llOQsl : bladman etahl llQod . I2OO1I : kllen et al.l I2OO2I : 
Petit et al. 2006 : Fraser et al. 20081 ). This suggests that the size distribution of KBOs is a power- 
law with logarithmic slope q ~ 4.5. The deepest survey of KBOs on the ecliptic, presented by 
Bernstein et al.l (|2004l ) , found a derth of faint objects demonstratin g that the LF must "roll-over" 
or break to a shallower slope at m(R) < 27 (jBernstein et al.ll2004l ). Interpreting the roll-over in 
terms of the Kuiper belt size distribution implies that the large-object power-law size distribution 
breaks at a KBO diameter ~ 100 km assuming 6% albedos. Models of ac cretion in the outer 
solar system predic ted such a break at least an order of magnitude smaller (jKenyon &: Bromley 



2001 



Kenvonll2002l ). 



The goal of our work was to confirm the results of iBernstein et al.l (|2004l ) , and accurately 
measure the shape of the LF at and beyond the roll-over magnitude. From this, the shape of the 
under-lying KBO size distribution can be inferred. We present here a survey of the ecliptic with 
limiting magnitude m{R) ~ 27, and use this along with previous observations to determine the 
shape of the LF. In section 2 we present our observations and in sections 3 we present our survey 
results. In sections 4 and 5 we consider past ecliptic Kuiper belt surveys and consider various LF 
functional forms. In sections 6 and 7 we present our analysis and a discussion of the results. 
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2. Observations, image processing, and characterization 

New observ ations were made in service mode with Suprime-Cam on the Subaru 8.2 m telescope 
(jMivazaki et al ]Eo02). Suprime-cam is a 10 chip mosaic CCD camera with a 34'x27' field-of-view 
(FOV) and a 0".202 pixel scale. Observations were made in the r-band on the nights of April 
22nd (night 1), and May 8th (night 2) 2008. The same target field was observed each night, 
targeting the ecliptic with coordinates a = 13'^46™57*.7, 5 = — 10°44 00 . Additionally, each night, 
a photom etric reference ima ge was taken of the Canada-France-Hawaii Supernova Legacy Survey 
D2 field (lAstier et al.M2006l ). Details of the observations are presented in Table [H Presented 
in Figure [1] is a summary of the image quality of both nights; as can be seen, the nights were 
photometric. 

Here we quickly describe our image processing techniques and subsequent moving object search 
for these data. A m ore thorough discussion of these techniques can be found in our previous work 
(jFraser et al.l 120081 ). We describe in more detail here the additional image processing steps that 
were necessitated by the use of Suprime-cam. 

All images went through the same image pre-processing before the image search step. From the 
overscan strip, the bias level was measured for each image, and subtracted. It was found that the 
bias level varied by ~ 100 ADU (10% of the typical bias value) along the bias strip. Proper removal 
of this variation required that a fourth order single-dimensional polynomial be fit in a least-squares 
sense to the overscan region collapsed to a single column, and was done on an image-by-image 
basis. This polynomial bias-fit was removed from the raw images. 

An average fiat-field was produced for each chip from the sky-flats provided from the service 
mode observer (10 from night 1, 6 from night 2) using the mask_mkflat_HA routine from the sdfred 



(Yagi et al. 


2002: 


Ouchi et al. 


2004) 



sky-flat images. Regions near detected sources are ignored when the average flat is produced. 

The Suprime-cam held is vignetted, with ~ 40% throughput-loss near field corners. This 
vignetting needed to be accounted for when adding artificial sources to the data as the artificial 
sources must accurately refiect the sensitivity variations across each chip. Thus, the average flat- 
fields for each chip were separated into a quantum-efficiency variation map, and an illumination 
pattern. The illumination pattern was created by running a 20 by 20 pixel median filter over the 
average fiat-fields. The quantum-efficiency variation map was set as the normalized difference of 
the flat-field and the illumination pattern, and was divided out of each of the over-scan subtracted 
images producing smooth images that revealed the vignetting pattern. 

The Suprime-cam FOV is largely distorted such that a Kuiper belt object (KBO), with a rate 
of motion > 3"/hr., whose motion would be linear on the sky, would not be linear in the images. 
Thus, s patial distort ions needed to be removed. Image distortions were measured with the dofit 
routine (iGwynl l2008l ) . This routine compares a reference source list containing true astrometric 
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positions of sources, and positions of those sources within the images; from this, it calculates a 
3rd order spatial distortion map. The image of the D2 field used to photometrically calibrate our 
images also provided a perfect astrometric reference, as sources in the field have been astrometrically 
calibrated with residuals better than 0.3". The distortion map was measured from the D2 field 
image of each ni ght, and then used to resample all images onto a spatially flat image using the 



SWarp package (jBertin et al.l 12002 ). Flux was conserved during this resampling, and produced 
images with a 0.1986" pixel-scale. The spatial rms residuals of the resampling were 0.06", or 1/3 
of a pixel as confirmed using the astrometric positions of stars in the USNO catalogs. 

The point-source image shape varied significantly across each chip, and warranted the use of a 
spatially variable point-spread-function (PSF). It was found that a PSF whose shape varied linearly 
with image position could accurately re-produce the point-source image shape variations. For each 
chip, a gaussian PSF with look-up table was generated from 15 bright, unsaturated, hand-selected 
point-sources. 

For each chip, a random number of artificial moving sources (between 150 and 250) were 
generated with random rates and angles of motion consistent with objects on circular Keplerian 
orbits between 25 and 200 AU, and with random fiuxes consistent with that of point sources between 
23 and 28.5 mag. These artificial sources were implanted blindly in the data, and the source list was 
revealed only after the search was complete. Additionally, 10 23rd mag. sources were implanted. 
These sources had flux sufficient to flag errors in the image combining algorithms. 

Artiflcial sources with sky motions larger than 0.2 pixels were split up into a number of dimmer 
sources with total flux equal to the original source. The centres of the dimmer sources were shifted 
to account for the motion of the object in the images. The flux of each artiflcial source was varied 
from image to image to match the average brightness variations of 20 reference stars with respect 
to a reference image (that with the lowest airmass). Additionally, the fluxes of the artiflcial sources 
were scaled to account for the vignetting apparent in each chip. This was done by implanting 
the artiflcial source into a blank image. The blank image pixel values were then multiplied by the 
illumination pattern. This was subsequently summed with the vignetted sky images. After artiflcial 
object implanting, the illumination pattern was removed from the images. The result was images 
with spatially smooth backgrounds containing artiflcial sources with spatially varying sensitivity 
that matched that of the images. 

No attempt was made to subtract a bias pattern from the raw images as no bias frames were 
provided from the service observations. The sdfred documentation however, claims that this step 



is uri necessary, as the bias level variations of Suprime-cam are low (jYagi et al.l l2002l : iQuchi et al 



200^). 



To maximize the searchable area of the images, stationary objects were removed by subtracting 
an image template from the vignetting corrected images. The image template for each chip was 
created from an average of the search field images usin g an artificial skepticism routine which places 



little weight on pixels far away from the pixel average (IStetsonlll989l ) . Every fifth image was used to 
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create the image template, reducing the presence of subtraction residuals created behind all moving 
sources, while ensuring a high quality template such that all statio nary sources w ere sufficiently 



removed. The image subtraction was done using psfmatch3 routine (jPritchetl l2005l ) 



Before final image stacks were produced, a high pixel mask was applied to the subtracted 
images to mask out any spurious hot-pixels and cosmic-ray spikes. The masked-subtracted images 
were spatially shifted to account for the motions of moving sources. A grid of shift rates (0.4 — 4.5 
"/hr.) and angles (±15" off the ecliptic) were considered which covered the range of motions 
consistent with bound solar system objects in the Kuiper belt on prograde orbits. In our past 
searches, we found that the search depth is not very sensitive to the choice in grid spacing between 
angles, but is very sensitive to the choice in spacing between rates. We chose a rate spacing small 
enough such that a moving source would never exhibit a trail longer than twice its seeing disk at 
the rate and angle which best approximated its apparent motion. Our grid consisted of 19 separate 
rates and 5 separate angles. 

The data from each night were photometrically calibrated to the SDSS wavebands b y com- 



parin g source flux measurements in the D2 held to those quoted in the Mega-pipe project (jGwvn 



20081 ). A source's magnitude in the Subaru data is given by rsub = —2.5 log (j) + Z where b is the 
source brightness in ADU, t is the exposure time, and Z is the telescope zeropoint. For each chip, 
a linear relation of the form rsub = C{g' — r')Mega — f'Mega used to transformation the r'-band 
magnitudes in the Mega-pipe project to Subaru r-band magnitudes. The colour term, C, and the 
zeropoint Z were fit in a least-squares sense to the source-fiux measurements in the data, and the 
Mega-pipe magnitudes using at least 40 sources on each chip. This procedure determined the best 
fit zeropoint to an accuracy of ~ 0.04 mag. The results are presented in Table [2j We found that 
chip 00 exhibits a 0.4 magnitude decrease in sensitivity compared to the o ther 9 chips, which is 



consi stent with other photometrically calibrated Suprime-Cam observations (jYoshida Nakamura 



20071 ). 



The transformations of the r' and g'-band Mega-pipe magnitudes to the r', and g'-bands 
in the SDSS filter set are given by r'j^i^g^ = r'soss - 0.024(5' - r')sDSS and 5^^^^ = g'soss " 
0.153(5' — r')sDSS- Using the colour term Cav = —0.048 which is the average of colour terms from 
Table [21 we find that the transformation between r-band Subaru and SDSS magnit udes is given 



by r^ pyy = rsub + 0.018(5' - r')sDSS- For typical KBOs (5' - r')sDSS ~ 0.7 mag (IFraser et al 



20081 ). Thus, the r-band Subaru and SDSS magnitudes differ by ~ 0.01 mag. which is smaller 
than the uncertainty of the telescope zeropoint. Hence, we use the approximation that for KBOs, 

^SD.SS — ''"Sub- 



3. Survey and results 



From the first night, the grid of image stacks for each chip was manually searched by one 
operator producing a candidate source list. Moving sources were identified by their appearance 
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in the images; a moving source is round (or nearly) in the image stack that best compensated 
for the source's apparent sky motion, and exhibits a ch aracteristic t r ail in other stacks, helping 
the operator distinguish moving sources from noise. See iFraser et alj (|2008l ) for a more complete 
discussion of this technique. 

Implanted artificial sources were identified in the candidate source list if the candidate source's 
location was within a few pixels of the artificial source's true centre. The detected artificial sources 
allowed us to measure our detection efficiency as a function of magnitude, vi^'sDSs)- represent 
the detection efficiency of the full field by the functional form 



vir'sDSs) 



1 - tanh -^MSS_ 



a 



(1) 



where r/max is the maximum efficiency, is the half-maximum detection efficiency magnitude, and 
g is approximately half the width of the fall from maximum detection efficiency to zero. This 
function was fit, in a least-squares sense, to the detection efficiency of artificial sources. The best 
fit parameters were {T]max,n,9) = (0.932 ± 0.008,26.86 ± 0.02,0.44 ± 0.03). The best-fit curve 
represents the measured detection efficiency well (see Figure [2]). 



We fit the measured efficiency with r/' 



' SDSS 



1 - tanh IsDssz 



1 _ tanh ISDS£_ 



which has been found to repre sent the detection efficiency of previous KBO surveys (jPetit et al 



20061 : iFuentes &: HolmanI |2008| ) . The best-fit curve is presented in Figure [2j As can be see, there 
is very little difference between the more complicated efficiency representation and that of Equa- 
tion [H and both describe the measurements equally well. Thus, we find that the more complicated 
function is not warranted or necessary. 

Flux measurements were made for all identified moving sources on the shifted image stack 
which contained the roundest image for that source. Magnitudes were measured in a 3-4.5 pixel 
radius aperture, and corresponding aperture corrections were determined from the image profiles 
of ten 23rd mag. artificial sources on each chip. By comparing the flux measurements of the 
detected artificial sources to their true values, we were able to characterize our fiux measurement 
uncertainties. For background limited sources, the uncertainty in a source's measured magnitude, 

7 — Z 

r, is given by Ar = 7 IO 2.5 where Z is the telescope zeropoint, and 7 depends on the telescope 
and the observations (|Fraser et al.l l2008l ) . We found that 7 = 1.33 best described the observed 
uncertainties. 

38 sources were not associated with artificial sources. 36 of these are identified as newly 
discovered KBOs. Two of these were identified as coincidental overlap of multiple poorly subtracted 
saturated galaxies. Both of these false detections have fiuxes below the 50% detection efficiency. 
This is consistent wi th our previous sur vey - a non-zero false candidate rate occurred for detection 
efficiencies t] < 50% (jFraser et al.ll2008l ). We therefore truncate our detection efficiency at the 50% 
threshold and ignore the 8 sources faint-ward of this level. 
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From the positions of the objects on the first night, orbits were calculated using fit-radec 



([Bernstein Sz Khushalanill200d ) and positions of the objects on night 2 were predicted. While not 
all of our objects had follow-up detections, we had a 100% follow-up rate for those 7 objects whose 
predicted positions fell in the FOV of the second night's observations. This is consist ent with the 



100% follow-up rate above the 50% detection threshold we report from our past survey (IFraser et al 



20081 ). We are thus confident that all detected non-artificial sources with magnitudes above the 



50% threshold are real KBOs. 

For each detected KBO, three image stacks were created (6 for those with follow-up) of the 
first half, middle half, and last half of the vignetting removed, non-subtracted images from a night's 
observations, using the rate and angle of motion that produced the roundest source image. From 
these, flux measurements were made. The final source magnitudes we report are the average of all 
the flux measurements made for that source off the image stacks where the source was at least a 
few FWHM away from nearby bright stars and galaxies. The results of our survey are presented 
in Table El 



4. Luminosity function 

Previous surveys have found that the differential LF (number of KBOs with magnitude m to 
m + dm per square degree) of bright KBOs {m{R) < 26) mag. is well represented by 



S(m) = ln(10)alO' 



a{m—mo ) 



(2) 



where S(m) is the number of objects with magnitude m (usually R-band) per square degree, a 
is the power-law "slope", and nio is the magnitude for w hich the sky density of objects with 



magn it udes < m is one p er square degree (see for example iTrujillo et al.l (j200ll ) ; iGladman et al 
(j200ll ): iPetit et al.l (120061')). In our previo us work, we found a = 0.65 and mo{R) = 23.43 for a 



broad range of surveys (jFraser et al.ll2008l ). iBernstein et al.l (j2004l ) presented the faintest Kuiper 
belt survey to date. Using a subs tantial amount of Hub ble observations 100 orbits) they achieved 
a search depth of m{R) ~ 28.5 (IBernstein et al.l 120061 ) and discovered three KBOs. The dearth of 



detections was at least a factor of ~ 10 below the number expected, requiring that the LF has a 
steep logarithmic slope for bright objects, that "rolls-over" to shallower slopes for fainter objects. 



Fraser et al.l (120081 ) presented a survey in which they searched 3 square degrees of the ecliptic 
to a depth of m{R) ~ 25.4 mag. and substantially increased the number of known detections 
suitable for a determination of the LF faint-ward of m{R) ~ 23.5. They found that the LF was well 
re presented by a power-la w given by Equation [2j They concluded that the lack of objects observed 
by IBernstein et al.l ([200J) must be caused by a sharp break rather than a broad roll-over, and that 
this could not occur for m{R) < 24.4. 



Additionally, IBernstein et al.l (|2004l ) examined the LFs of the so-called cold and excited KBO 
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populations, i < 5° and i > 5° respectively. They found that the cold population exhibited 
a steeper bright object slope than the excited population. The statistical support for such a 
difference however, was not strong. Simulations that test different Kuiper belt formation scenarios 
suggest that a difference in the SDs of different KBO populations is caused by dif ferent formation 



condit ions /histories for those objects. Thus, a confirmation of the findings by iBernstein et al 



(|2004l ) is important, as it provides a very strong constraint on the Kuiper belt formation history. 



The survey reported here detected 28 KBOs that can be used to characterize the luminosity 
function (LF) and SD of KBOs. With these detections, we wish to address three questions: 

1. At what magnitude is the large object LF no longer a power-law? 

2. Is the observed break consistent with the LF expected if the SD has a sharp break at some 
object radius? If not, what is the true shape of the LF? 

3. Is the SD of the excited population different from that of the cold population? 



4.1. Luminosity Function forms 



In determining the correct LF shape, we consider three functional representations of the dif- 
ferential LF. The first is the power -law given by Equation El The second is the rolling power-law 
suggested bv Bernstein et al.l ((200^), given by 



E(m) = 5]23l0"('"-23)+-'(™-23)' 



(3) 



where a is the bright object slope, a' is the derivative of the logarithmic slope, and E23 is the 
number of objects per square degree at 23rd mag. 



The third functional form we consider was first presented by iFraser &: Kavelaara ([2003) • They 
showed that, if the SD is a power-law with slope for large objects that has a sharp break at 
objects with diameter Db to slope q2 for small objects, then the cumulative LF has the form 



N {<m) = A< 



ailO"!'" if m < rub 

6il0°2™-^62l0"i'" + 63(m) if mb<m< 1712 
ci + C2lO"2™ if m> 1712. 



(4) 



where S(m) 



dN{<m.) 

dm ' '^1 



gi-l 



and 02 



92-1 



In Equation m mfe = K + 5 log ^) (break 



5 """^ - 5 

magnitude) and m2 = K + 5 log (r|D^^) are the magnitudes for which objects smaller than the 
break diameter are detectable by a given survey for heliocentric distances larger than the Kuiper 
belt inner edge ri but smaller than the outer edge r2- 
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The parameter is a constant relating the diameter of an object to its apparent magnitude. 
Choosing a value of K is equivalent to choosing a constant albedo for KBOs. It is apparent that 
the lar gest objects have albed os as high as ~ 90 %, with smaller objects exhibiting much lower 
values (jStansberrv et al.l 120071 ) . Clearly the assumption that albedos are constant is incorrect, 
at least when considering the largest KBOs. Insufficient knowledge of the albedos for smaller 
D < 500 km KBOs - the siz e range which we probe in this work - exists for this assumption to 
be tested. Fraser et al. (j2008l ) discuss the effects of a varying albedo when a constant is assumed, 
and demonstrate that the inferred size distribution slope is incorrect. This effect however, is small 
compared to the uncertainties typical in the LF determination, given the available data. 

It is apparent, at least for the observe d sources, that an albedo of 6% is a typical value for 
D < 500 km KBOs (IStansberry et al.l 120071 ). and is the value we assume here. This sets K ~ 18.4 
as the R-band magnitude of a D ^ 1600 km object at 40 AU with a 6% albedo. This apparently 
typical albedo might very well be an observational bias. Until more data are available, it will not 
be known whether 6%, or a much larger value similar to large objects, is a more appropriate choice. 
Thus, any conclusions about the inferred break size come with this caveat; the real break diameter 
could be quite different dependening on the true albedos. 

This functional form has the advantage that the LF can be interpreted in terms of LF parame- 
ters («!, 02, rrifc, log A) or SD parameters {qi,q2, D^, log A). Both provide equivalent representations 
of the LF shape, and, assuming an albedo, the SD is directly inferred from the best-fit LF. 

The coefficient s a,;, b j, and q of Equation H] are functions of the KBO radial distribution. 
Fraser fc Kavelaara (|2008l ) have shown that the LF shape for < m < m2 is strongly dependent 
on the assumed KBO radial distribution, but the SD slopes qi and ^2, and the break diameter 
however, can be accurately inferred if data exists sufficiently far away from the break magnitude. 
They assume that the radial distribution is given by i VfA) (x A~'^, and sugge st that c = 10 is 
a good approximation to the radial fall-off observed by iTrujillo &: Brownl (|200ll ) . We assume the 
same radial distribution here. 



These three functions provide a good sample of LF forms to consider as they progress in 
levels of complication by a single degree of freedom - the parameters for the different LF forms are 
(a, mo) for Equation^ (Q;,a',Il23) for Equation^ and {ai,a2,miy, A) for Equation^ Thus, it is 
a straight-forward process to determine the number of parameters necessary to describe the true 
shape of the LF using simple statistics techniques. 



4.2. Data-sets 



In our previous estimates of the LF, we considered observations from surveys f or which the 
detect ion efficiency was well characterized as a function of magnitude. The F08 sample 



2001 



(120081 ) consists of n e ar ecl ip tic KBO su r veys p resent ed in 



); lAllen et al.l (120021 ): IPetit et al.l (|2006l ). and iFraser et all (120081 ) along with the data from 



Jewitt et al. 



Fraser et al 



(119981 ) : lGladman et al.l (|l998l . 
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Truiillo et alj (|200ll ) which we subdivided by detection efficiency and astrometri c position. The 



F08b sampl e is the F08 sample, plus the inclusion of the survey presented by iBernstein et al 
(j2004l . 120061 ) and the survey presented here. 



To measure any differences in the LF between the dynamically excited and cold populations, 
we define the F08bj<5 and F08bi>5 samples as those objects with inclinations smaller than and 
greater than 5° respectively. For these samples, we only consider surveys in which the majority of 
objects have been observed on arcs longer than 24 hrs, inclinations measured from the observations 
spanning only one night are not accurate enough to make an inclination cut reliable, and a large 
contamination between the two samples would occur . Thu s we exclude from the F08bi<5 and 
F08bi>5 subsamp les the Uran u s field from lPetit etaP (|20od l. the N11033, N10032W3, UNE, and 
UNW fields from lFraser et al.l (|2008l ). and the survey presented here. 



While an inclination cut does not directly probe the different KBO p opulations, ie. cold/hot- 
classical objects, resonators, scattered members, etc. (jGladman et al.ll2008l ). it does provide a crude 
means of separating those KBOs with 'highly' excited orbits from those with les s excited orbits 
The division of i = 5° is chosen to ease comparison with the LF presented by IBernstein et al 



(|200J). For the F08bj<5 and F08bj>5 samples, no heliocentric distance cut was made. 



Note: the analysis f or this manuscript was alr eady complete before the authors were aware of 
the survey presented by lFuentes &: HolmanI (120081 ). The moderate increase in the total number of 
detections, ^ 20%, will not substantial ly alter the results presented here. Thus, we did not include 
the survey of lFuentes &: HolmanI (j2008l ) in our analysis. Their results however, are discussed below. 



4.3. Analysis 



In determining the corr ect LF shape and be s t-fit parameters of the vario us functional forms, we 
adopt the same practices as lFraser et al.l (|2008l ): iFraser &: Kavelaara ()2008l ). These are as follows: 



1. We cull from each survey all sources faintward of the 50% detection efficiency of that survey, 
and set the detection efficiency to zero faintward of that point. 



2. We of fset all magnitudes to R-band using the average KBO colours presented in lFraser et al, 



3. We fit the differential LF, S(m) 



dN{<m) 
dm 



using a maximum likelihood technique. 



4. We adopted the same form of the likelihood equation as presented in iFraser et al.l (120081 ) 
which treats calibration errors and variable sky densities as separate nuisance parameters. We 
marginalized the likelihood over these parameters before determining the maximum likelihood. 



5. We adopted the s ame prior ranges of the colour and density parameters as presented in 
Fraser et al.l (|2008l ) (±0.2 magnitudes for the colour parameters, and ±0.4 for the logarithm 
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of the density parameters) . 
In determining the quahty of the fits, we utihze the Anderson-Darhng statistic 

Jo P{m) (1 - Pirn)) * ' ' ' 

where P{m) is the cumulative probabihty of detecting an object with magnitude < m and S{m) 
is the cumulative distribution of detections. We calculate the probability, P(A > Aq^s) of finding 
a value of the Anderson-Darling statistic. A, larger than that of the observations for a given LF 
function and parameter set by bootstrapping the statistic, is. randomly drawing a subsample of 
objects from that LF equal to that detected in a particular sample, and fitting the LF to that 
random sample and computing A. In the random sampling process, we include random sky-density 
and colour offsets which represent those expected from the real observations. Values of P(A > Aohs) 
near indicate that the functional form is a poor representation of the data. 

To determine if the more complicated LF functional forms are statistically warranted, we utilize 
the log-likelihood test = ~21og-2^. The logarithm of the ratio of maximum likelihoods of the 
simple and more complicated functions, L and L', are distributed as a chi-squared variable with a 
number of degrees of freedom equal to the difference in free parameters between the two functional 
forms. A table of values can be used to determine the signifi cance of the improvement in the 



LF fit between the two functional forms. (jKotz Johnsonlll983l ). 



5. Results 
5.1. F08b Sample 

Presented in Table [Hare the best-fit parameters for the three different LF forms, and P(A > 
^obs) of each fit. The best-fit power-law slope to the F08b data- set, a = 0.5 8 , is m oderately 



shallower than, but is consistent with the best-fit a = 0.65 from iFraser et al.l (|2008l ). The fit 



however, is a poor description of the data as evidenced by the low P(A > Aobs) < 0.02. 

The best-fit of the rolling power-law to the F08b data-set is (a,a',S23) = (0.8,-0.06,0.82). 
The maximum likelihood value of this fit is increased by more than two orders of magnitude over 
the power-law which warrants the inclusion of a third degree of freedom, and the fit is an acceptable 
description of the F08b sample, with P(A > Aof^g) = 0.4. 

The best-fit of the LF given by Equation H {log A, ai, 02, rrib) = (23.56,0.76,0.18,24.9), is an 
acceptable fit to the data with P(A > Aobs) = 0.4. The maximum likelihood value is increased by 
more than two orders of magnitude over the power-law with a log-likelihood chi-square, = 11.7, 
and by a factor of 2 over the rolling power-law with = 1-4. We find that the best-fit from 
Equation m is preferred over the power-law at greater than the 3-sigma level. This fit is preferred 
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over the rolling power-law at ~ 80% significance. Thus we find that both the rolling power-law 
and broken power-law of Equation [J] provide equally adequate descriptions of the F08b sample. 
Note: the Monte-Carlo simulations we have done to calculate P(A > Aobs) have also demonstrated 
that on average, the best-fit LF parameters determined from our maximum likelihood technique 
reproduce the input parameters of the simulations for all LF functionals forms we have considered 
here. 



When the observations of lBernstein et al.l (|2004l ) are excluded from the FOSb sample, the best- 
fit power-law has para meters {a, mo) = (0.65,23.42), which is nearly identical to that found from 
bv lFraser et al.l (|2008l 1. The power-law is a moderately adequate description of the observations 
with P(A > Aofjg) =0.1. This result hints that the Kuiper belt LF exhibits a break within the 
magnitude range of the FOSb sample, but fits to the survey data with m{R) < 27 do not explicitly 
require a broken or rolling power-law description. 

Presented in Figures [3l [H and [5] are the credible regions of the fits to the FOSb sample. 
The function given by Equation [J] is derived from a size-distribution similar to that expected from 
accretion calculations, and thus provides a proper interpretation of the under-lying size-distribution. 
We find that the LF is consistent with a broken power-law size-distribution with large and small 
object slopes = 4.8 ± 3 and ^2 = 2 ± 2, and a break diameter, = 62 it 40 km assuming 6% 
albedos. 

Presented in Figure[6]is a histogram of the FO Sb sample, the best -fit LFs for the three functional 
for ms we consider, the best-fi t power-law LF from lFraser et al.l (|2008l ). and the best-fit LF presented 
bv lFuentes k. HolmanI ((20081). As can be seen, power-laws do not provide acceptable fits to the data 
with m{R) ~ 28, and a more complicated function is required. Equations [3] and H] both provide 
reasonable descriptions of the observations. 

The user is cautioned from drawing any further conclusions from Figure [H The FOSb sample 
contains data from many different surveys. Thus the sample includes any calibration errors and 
sky-density variations be tween survey s whic h might lead the reader to a false conclusion about the 
structure of the LF (see iFraser et al.l (|2008l ) for a discussion of the magnitude and significance of 
these effects). 



The best-fit LF presented by lFuentes &: HolmanI (|2008l ) is a harmonic mean of two power-laws 



with bright-end slope ai = 0.7 and faint-end slope 02 = 0.3 with a break magnitude of Req ~ 24.3. 
These parameters are consistent with the slopes and break magnitude of our best-fit non-power-law 
LFs. They find however, a factor of 3 higher sky-density which is inconsistent at greater than the 
3-sigma level with the range of sky densities we find for similar bright-end slopes (see Figure H]) . 
The source of the increased number of detections compared to that expected from fits to the FOSb 
sample is unknown. 
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5.2. F08b,;>5 Sample 

When we consider the F08bi>5 sample, we find that the best-fit power-law has slope a = 0.55 
with P(A > Aobs) < 0.04, and is a poor description of the data. The sample has no objects fainter 
than m{R) ~ 25.8. The lack of faint objects requires a break near this point. The best-fit rolling 
power-law with parameters (a,a',S23) = (0.74,-0.06,0.46) is an acceptable fit to the data with 
P(A > ^obs) — 0.1. The best-fit broken power-law model is an acceptable description of the data 
with P(A > lS.oBs) = 0.1. The log-likelihood ratio implies that the broken power-law model of 
Equation S] is preferred over the rolling power-law at the 90% significance level. The best-fit has 
ai = 0.7 and breaks at D[, = 36 km. Because the F08bj>5 sample has no detections faint-ward of 
the break, both the break diameter and small-object slope of this sample are poorly constrained. 
The maximum likelihood is found with a break at m{R) = 25.8 with no objects faint-ward of this 
point, ie. a2 = —oo. Note: we consider slopes, 02 > —6 to avoid numerical errors in the likelihood 
calculations. 

Presented in Figures [71 EJ and [9] are the likelihood contours of the best-fit of Equation [4] to the 
F08bj>5 sample and the histogram of those data. As can be seen, the LF is well described by a 
power-law for m{R) < 25.8. The lack of faint objects requires that a break must exist at D^, < 100 
km. 

The faint-object slope is highly uncertain, with a 1-sigma upper-limit of 02 ~ 0.6. This fit 
however is consistent with that of the F08b sample best-fit at the ~ 1-sigma level. Though the 
inclinations of the detections we present in this survey are inaccurate, they hint that the break in 
the LF for the F08bj>5 sample is not as sharp as the fit would suggest. The detections from this 
survey are consistent with the same break diameter (36 km) that breaks to a much flatter faint 
object slope similar to the 1-sigma upper-limit of the fit (02 ~ 0). The fit to the F08bj>5 sample 
however, constrains the break to be brightward of m[R) < 26.8. Additional data brightward of 
m{R) ~ 27.5 with accurate inclinations is required before the break location can be accurately 
constrained. 



5.3. F08bi<5 Sample 

The low-inclination LF is sufficiently described by a power-law. The best-fit power-law is an 
acceptable description of the sample, with parameters (a, nio) = (0.59, 24.0) and an Anderson- 
Darling statistic, P(A > Aobs) = 0.15. Indeed the best-fit of Equation H] is found when the break 
occurs at the bright-end limit of the data, suggesting there is no strong evidence for a break in the 
F08bi<5 data. 

Presented in Figures fTOl and [TT] are the likelihood contours of the fit to the low-inclination LF, 
and the histogram of the observations. As can be seen, there is no apparent evidence for a break 
in the magnitude range of the observations. 
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To test whether or not the excited and cold samples have different LFs, we utilized the 
Anderson-Darling test described above (see Section I4.3P to determine whether or not the bro- 
ken LF fit of the F08bj>5 sample is an adequate description of the F08bj<5 sample. We however, 
set a2 = rather than the best-fit value, as this is a more realistic faint object slope than that of the 
best-fit of the F08bj>5 sample {02 ~ —00), as suggested by our survey. This slope is still within the 
1-sigma credible region of the fit (see Section [5.2p . We found that for the F08bj<5 sample, an LF 
with {ai, a2,mb,log A) = (0.7,0,26.2,22.7) has P(A > Aobs) = 0.4, and is a sufficient description 
of the cold sample LF. We therefore conclude that we have no evidence that those KBOs with i > 5 
have a different LF than those with i < 5. 



Bernstein et al.l (120041 ) found that the LFs of the cold and excited populations were different. 
Fuentes fc HolmanI (j2008l ) found similar results for their survey alone, and concluded that the low 
inclination population LF exhibits a bright-end slope a ~ 0.8 — 1.5, steeper than ai ~ 0.7 found 
when they considered objects of all i in their survey. When they considered all available data - the 
F08b sample, but not including the survey presented here, and including their own search - they 
found that the low-inclination group exhibited a significantly steeper bright-object slope than the 
high-inclination sample. This finding is intriguing, as we find no evidence for a difference in the 
SDs of the low and high inclination groups. It is possible that their results are a consequence of 
not considering calibration and density offsets when performing the likelihood calculations, as we 
have done here. Clearly however, more observations are needed which probe the entire range of 
current observations before these results can be clarified. 



6. Discussion 



The best-fit parameters for the LF presented in Equation [3] imply that the size distribution is 
a power-law with slope qi ~ 4.8 for large objects, which breaks to a shallower slope at diameter 
Di, ^ 50 — 95 km assuming 6% albedos. Comparison of this size distribution to models which evolve 
a population of planetesimals and track their size distribution as a function of time in the Kuiper 
belt re gion can place a constraint on the forma t ion processes and the dur ation of accretion in that 
region (jKenvon &: Bromlevll200ll : lKenvonll2002l : iKenvon Sz Bromlevll2004l ). 



In the early stages of formatio n, run-a,wav g rowth occurs, and very rapidly grow objects as large 
as 10^ km in the Kuiper belt ( Kenyonll2002l ). During this process, a steep-sloped large object 
size distribut ion deve l ops, w hich flattens with time as more objects become "large". Calculations 
presented in iKenvonI (j2002l ). which simulate planet accretion for conditions in the Kuiper belt 
expected for the early solar system, imply that, in the absence of influences from Neptune, the 
modern-day large object slope would be shallower than that observed if accretion lasted the age 
of the solar system. They find that for an initia,! Kuip er belt mass similar to that predicted from 
the minimum-mass solar nebula model (|Havashilll98lh - much more massive than the current belt 
(jFuentes fc HolmanI l2008l ) - accretion must have been halted after 100 Myr. If KBOs are weak 
(easier to disrupt), then accretion might have gone on for as long as 1 Gyr (see Figures 9, and 10 
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from lKenvoiill2002l ). Some event (s) must have halted this process by clearing the majority of initial 
mass out of the belt before the slope became too shallow. 



Kenvon &: BromlevI (|2004l ) has demonstrated that weaker bodies will exhibit a size distribution 



with a larger break diameter Db than stronger bodies would if they underwent the same evolutionary 
history. They calculated the size distribution expected from a Kuiper belt under the gravitational 
influence of Neptune which evolved for the age of the solar system. In this model, break diameters 
as large as 60 km were produced only for the weakest bodies they considered. This model however, 
produced a large object slope much too shallow to be consistent with the observations, implying 
that accretion over the age of the solar system cannot have occurred. 

The existence of such a large break diameter in the KBO SD implies that, KBOs must be quite 
weak (strengthless rubble piles), and have undergone a period of increased collisional evolution. This 
suggests a scenario in which the event responsible for the early end to accretion, and the clearing 
of the majority of the mass in the Kuiper belt, also increased the rate of collisional evolution for a 
time, pushing the break diameter to the large value we see today. 

Collisional disruption would be subst antially increased for a p eriod of time during the scat- 
terin g of planetesimals by a rogue plane t (jGladman &: ChanI l2006l ) . or during a close stellar pas- 
sage (llda et al.ll200Cll : iLevison et al.ll2004l ). During these scenarios, relative velocities of KBOs are 
increased, causing collisions between bodies to result in catastrophic disruption, rather than accre- 
tion as would otherwise be the case. The combination of scattering and collisional grinding would 
produce a substantial mass loss in the Kuiper belt region, and produce the large knee observed 
currently. Increased bombardment would also occur if Neptune migrated outwards. In such a sce- 
nario, Neptune formed closer to the Sun than its current location. Gravitational scattering of small 
planetesimals transfer angular momentum to Neptune, causing the gas-giant to migrate outwards. 
During this process, some planetesimals have their orbits excited, and are thrust into the modern 
day Kuiper belt creating the orbital distribution we see today. The remaining planetesimals are 
cleared from the system during close encounters with Neptune causing a rapid and substantial 
mass-loss during the migration process. Various incarnations of this attractive scenario can provide 
the necessary collisional bombardme nt of KBOs, and simultaneously account for some of the other 



Levison et al.ll2008l ^. 



observed features of the K uiper belt (jMalhotralll993l : lLevison &: Morbidelh 
2005 



2003 



Hahn fc Malhotra 



The Neptune migration scenario predicts that the more excited population will have a more 
evolved size distribution, with a shallower large-object slope - objects which originated closer to 
the Sun have stronger encounters with Neptune, and thus have the most excited final orbits. These 
now excited objects would have originated from a more dense region, and hence underwent more 
rapid accretionary evolution, than objects which originated further from the Sun producing the 
shallower large object slope. We do not see any differences in the large object slopes of the cold 
and excited populations. Our findings suggest that, assuming migration of Neptune occurred, the 
total distance travelled by Neptune was sufficiently small such that the rate of planet formation 



and hence the size distribution of all objects scattered by Neptune and implaced in the Kuiper belt 
could not be substantially different. 



7. Conclusions 

We have performed a survey of the Kuiper covering ~ 1/3 a square degree of the sky using 
Suprime-cam on the Subaru telescope, to a limi ting magnitude oi nif^oiR) ^ 26.8 and have found 



36 new KBOs. Using the likelihood technique of lFraser et al.l (|2008l ) which accounts for calibration 



errors and sky density variations between separate observations, we have combined the observations 
of this survey with previous observations, and have found that the luminosity function of the Kuiper 
belt is inconsistent with a power-law with slope ai = 0.75, but must break to a shallower slope 
a2 ~ 0.2 at magnitudes m{R) ~ 24.1 — 25.3. The luminosity function is consistent with a size 
distribution with large object slope qi ~ 4.8 that breaks to a shallower slope ^2 ~ 1-9 at a diameter 
of 50 — 95 km assuming 6% albedos. We have found no conclusive evidence that the size 
distribution of KBOs with i < 5 is different than that of those with i > 5. 
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Table 1. Observation Details. 



Night 



Date 



Seeing (" ) Number of Exposures 



a 



S 



2 



1 



April 22 2007 
May 8 2007 



0.4 - 0.7 55 X 200 s. 

0.4 - 0.8 53 X 200 s. 



13:46:57.7 -10:44:00 
13:46:57.7 -10:44:00 
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Table 2. Data Calibration Details. The zero-point, Z, and colour term, C, of each chip of 

Suprime-Cam (see Section 2). 



Chip 


Z (mag.) 


C 


00 


27.48 


-0.040 


01 


27.81 


-0.050 


02 


27.78 


-0.065 


03 


27.85 


-0.020 


04 


27.73 


-0.095 


05 


27.89 


-0.020 


06 


27.82 


-0.045 


07 


27.83 


-0.040 


08 


27.86 


-0.040 


09 


27.75 


-0.065 
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Table 3. Survey Results. rsDSS is the average r-band magnitude, n is the mimber of flux 
measurements used in the average, A is the object's hehocentric distance, and i is the object's 

inchnation. 



Object 


rsDSS (mag.) 


n 


A 




Follow- 


c9a2 


24.038 ± 0.0424 


3 


40 ±3 


22± 14 


n 


cla3 


24.368 ± 0.0575 


3 


44 ±3 


9± 12 


n 


c6a3 


24.394 ± 0.0585 


6 


45 ±2 


5±2 


y 


clal 


24.991 ± 0.1018 


3 


44 ±3 


1±7 


n 


c3al 


25.179 ±0.1212 


3 


42 ±3 


0±2 


n 


c5a2 


25.245 ±0.1293 


3 


41 ±3 


3.3 ±10 


n 


c5a5 


25.373 ±0.1444 


3 


49 ±3 


3± 14 


n 


c8a2 


25.397 ±0.1485 


3 


47 ±4 


15± 15 


n 


c4al 


25.432 ± 0.1527 


3 


45 ±3 


9± 13 


n 


c9al 


25.434 ± 0.1527 


3 


35 ±3 


28± 16 


n 


c9a4 


25.438 ± 0.1541 


3 


48 ±3 


2± 13 


n 


c8al 


25.458 ± 0.1569 


3 


43 ±3 


2± 11 


n 


c2al 


25.462 ±0.1569 


3 


48 ±3 


7± 14 


n 


c4a3 


25.569 ± 0.1737 


3 


38 ±3 


5± 10 


n 


c4a5 


25.823 ± 0.2187 


3 


42 ±3 


5± 11 


n 


c5a9 


25.875 ± 0.2311 


3 


45 ±3 


5± 12 


n 


c6al 


26.019 ± 0.2629 


6 


37 ±3 


2±9 


y 


c0a2 


26.068 ± 0.2753 


6 


45 ±2 


1.1 ±0.3 


y 


c2a3 


26.279 ± 0.3340 


4 


43 ±3 


5± 12 


y 


cla4 


26.282 ± 0.3340 


3 


44 ±3 


4± 12 


n 


cOal 


26.334 ± 0.3498 


6 


45 ±3 


14 ±5 


y 


c6a2 


26.382 ± 0.3663 


6 


44 ±2 


3±1 


y 


c2a2 


26.583 ± 0.4404 


6 


39 ±2 


10 ±4 


y 


c4a6 


26.610 ± 0.4527 


3 


43 ±4 


23 ±17 


n 


c4a2 


26.673 ± 0.4784 


3 


45 ±3 


12± 13 


n 


cla5 


26.680 ± 0.4828 


2 


37 ±3 


27± 17 


n 


c5a3 


26.688 ± 0.4873 


2 


44 ±3 


7± 12 


n 


c5a8 


26.737 ± 0.5103 


3 


47 ±3 


7 ±13 


n 


c5a7 


26.858 ± 0.5699 


3 


44 ±3 


1 ±9 


n 


c4a7 


26.859 ± 0.5699 


3 


34 ±3 


9±9 


n 


c9a3 


26.883 ± 0.5805 


1 


47 ±5 


75 ±29 


n 


c7al 


26.925 ± 0.6079 


3 


45 ±3 


16± 15 


n 


c5al 


27.120 ± 0.7241 


3 


43 ±3 


1±8 


n 


c5a6 


27.276 ± 0.8391 


3 


47 ±3 


6± 14 


n 


cla2 


27.353 ± 0.8950 


2 


55 ±4 


12± 18 


n 
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Table 3 — Continued 
c5a4 28.807 ± 3.4344 3 42 ±3 Oil n 



Table 4. Luminosity Function Fits (LF and SD parameters) and associated Anderson-Darling 
statistics. < indicate 1-sigma upper limits. We do not present fits utilizing functions with more 
degrees of freedom than are statistically warranted by the observations. The F08bi<5 sample was 
not fit with Equation [J] but is adequately described by LF parameters within the 1-sigma range of 
the fit of the F08bj>5 sample. NOTE: q2 = — oo represents a complete absence of objects smaller 

than the break diameter. 



Data-set Power-Law Rolling Power-law Broken Power-law 



(a, mo) P(A>Aob,) (S23,a,a') P(A > A^f,,) (log A, ai, 02, "ife) P(A > A, 

{\ogA,qi,q2,Db) 



F08b 


(0.58,23.31) 


< 0.04 


(0.82,0.80, 


-0.06) 


0.4 


(23.6,0.76,0.18,24.9) 
(23.6,4.8,1.9,62) 


0.4 


F08bi<5 


(0.59,24.0) 


0.15 












F08b,>5 
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Fig. 1. — Reference star magnitude versus airmass (top) and full-width-half-maximum versus ob- 
servation time (bottom) for the April 22nd (left) and May 8th (right) nights. Typical uncertainty 
in the the magnitudes is shown. The reference star displays a ~ 0.1 mag. variation in each night, 
which is smaller than the typical uncertainty in the flux measurements of each of our sources. 
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Fig. 2. — Detection efficiency histogram with 1-sigma poisson errorbars, and the best-fit efficiency 
functions {vmax,r^,g) = (0.932,26.86,0.44) (soUd) and {Vmax,r*, g, g') = (0.941,27.14,0.333,0.758) 
(dashed-dotted) . 
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Fig. 3. — The 1,2, and 3-sigma credible regions of the best-fit roUing power-law to the F08b sample 
with E23 held at its best-fit value. 
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Fig. 4. — The 1,2, and 3-sigma credible regions of the best-fit of EqlHto the F08b sample with 0:2 
and Df, held at their best-fit values. 




Fig. 5. — The 1,2, and 3-sigma credible regions of the best-fit of Eq|3]to the F08b sample with ai 
and A held at their best-fit values. NOTE: This figure assumes that KBOs have 6% albedos (see 
Section 4.1). 
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Fig. 6. — Histogram of the data co nsidered in the F08b sample with al magnitudes shifted to R- 
band using the colours presented in iFraser et al.l (|2008l ). Errorbars are 1-sigma poisson intervals. 
The solid curve is the best-fit LF of Eq uation [H T h e das hed curves are the best-fit power-law 
LFs from this manuscript (shallow) and iFraser et al.l (|2008l ) (steep). The das hed-dotted curve is 
the be st- fit rolling power-law. The dotted curve is the best-fit LF presented by lFuentes &: Holman 
(|2008l ). Object diameters (km) are given assuming 6% albedos at 35 AU. 
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Fig. 7. — The credible regions of the best-fit of Equation [4] to the F08bj>5 sample. Contours are 
the 1, 2, and 3-sigma credible regions of the fit with 02 and D\, held at their best-fit values.. 
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Fig. 8. — The credible regions of the best-fit of Equation [J] to the F08bj>5 sample. Contours are 
the 1, 2, and 3-sigma credible regions of the fit with ol\ and A held at their best-fit values. The 
dotted horizontal line marks the best-fit large object slope, and demonstrates the need for a break 
to shallower slopes for this sample. NOTE: This figure assumes that KBOs have 6% albedos (see 
Section 4.1). 
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Fig. 9. — Histogram of the data considered in the F08bj>5 sample. Errorbars are 1-sigma poisson 
intervals. The solid curve is the best-fit of Equation HI The dashed curve is the LF expected 
when ai, Df, and A are held to their best-fit values, and 02 is taken at the upper 1-sigma limit, 
and represents the l - sigma upper limit to the LF set by the absence of detections with i > 5 
from IBernstein et al.l (|2004l ) . The stars are the histogram of the Subaru survey presented in this 
manuscript. These points however, are unreliable as the inclinations for the objects discovered in 
this survey are highly uncertain. Object diameters (km) are given assuming 6% albedos at 35 AU. 
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Fig. 10. — The credible regions of the best-fit power-law to the F08bj<5 sample. Contours are the 
1, 2, and 3-sigma credible regions of the fit. The best-fit power-law is consistent with the coUisional 
equilibrium slope q = 3.5 and is a statistically sufficient description of this data-set. The best-fit 
broken power- law of the F08bj>5 sample is a sufficient description of this data-set. 
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Fig. 11. — Histogram of the data considered in the F08bj<5 sample. Errorbars are 1-sigma poisson 
intervals. The solid curve is the best-fit power-law and is an adequate description of the -F086j<5 
sample. Object diameters (km) are given assuming 6% albedos at 35 AU. The dashed line is a broken 
LF from Equation [4] with {ai,a2,'mb,logA) = (0.7,0,26.2,22.6) and is an adequate description of 
the data. The dashed curve is consistent at the 1-sigma level with the best-fit LF of the F08bi>5 
sample, demonstrating that, statistically speaking, the -F086j<5 and -F086j>5 samples have the same 
luminosity functions. 



